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Abstract 

A theoretical description of electromagnetic waves in the back- 
ground of a (weak) gravitational wave is presented. Explicit expres- 
sions are obtained for the Stokes parameters during the passage of a 
plane-fronted gravitational wave described by the Ehlers-Kundt met- 
ric. In particular, it is shown that the axis of the polarization ellipse 
oscillates, its ellipticity remaining constant. 
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1 Introduction 



In recent years there has been considerable interest in the detection of gravita- 
tional waves using interferometric methods pQ . Accordingly, it is convenient 
to have a precise description of how an electromagnetic (EM) wave is affected 
by the passage of a gravitational wave. The goal of the present article is to 
describe this interaction solving Maxwell's equations in the corresponding 
space-time background. For this purpose, the metric of Ehlers and Kundt 
[2] provides an appropriate description of a plane-fronted gravitational wave 
and will be used throughout the paper. Particular attention will be focused 
on the Stokes parameters, as these are directly observable and provide a 
complete description of the EM waves. 

In section 2 of this article, Maxwell's equations in the background of the 
Ehlers-Kundt metric are solved for a weak gravitational field; all relevant 
formulas for the EM field are then deduced within the short wave-length 
approximation. Section 3 is devoted to the calculation of the Stokes param- 
eters describing the polarization of the EM wave. The main physical result 
is that the eccentricity of the polarization ellipse remains constant but its 
axis rotates during the passage of a gravitational wave. If the latter has an 
oscillatory behavior, it will produce an oscillation of the polarization ellipse 
axis with the same frequency. 

2 Plane-fronted gravitational waves 

The metric of a plane-fronted gravitational is given by [2] 

ds 2 = g^dx^dx" = H(u, x, y) du 2 — 2du dv + dx 2 + dy 2 , (2.1) 

where 

H(u, x, y) = a(u)(x 2 - y 2 ) + 2b(u)xy , (2.2) 

and a and b are functions of the null-coordinate u. The inverse metric tensor 
is 
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(2.3) 



Here and in the following, coordinates are taken in the order (u, v, x, y). 
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Minkowski space-time in Cartesian coordinates (t,x,y,z) with signature 
(— , +, +, +) is recovered if H = 0. In this limit, one can identify the coordi- 
nates as 

1 , 1 / 

U=-j=(t-Z), v = -j=(t + z). 

As for the D'Alembertian operator, it takes the form 

d ( ..„ d 



□ 
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d 2 d 2 d 2 d 2 

du 8v i ■> iV) Q v 2 Q x 2 Qy2 ( ) 

(notice that y/^g — 1). 

The electromagnetic field tensor F a/3 can be obtained from a potential A a 
such that 

F a() = d a A fi - dpA,,, (2.5) 
and must satisfy the Maxwell equations in vacuum, 

° F Q/3 = 0. (2.6) 



dx? 



A convenient choice for the gauge condition is A v = 0, since in this gauge the 
components of decouple and the following explicit forms of the equations 
are obtained: 

OA x — — UA y , (2.7) 

UA U = d v (A x d x H + A y d y H), (2.8) 

as can be checked by direct substitution. The last of these equations is 
equivalent to 

d v A u - d x A x - d y A y = 0, (2.9) 

which is consistent with Maxwell's equations. 

The general solution of the equation OA = 0, valid to first order in the 
gravitational potentials, is of the form 

A oc exp{ik^ + (2.10) 

where 

$(u, x, y) = --k v [a"(x 2 - y 2 ) + 2(3"xy 
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- (a'k x + P'ky) x + (-/3'k x + a'ky) y-^- \a(k 2 x - k 2 y ) + 2(3k x k y , (2.11) 
with 



a{u) = a"'(u), b{u) = f3"'(u) , 

the primes denoting derivative with respect to u. Here and in the following, 
all quadratic or higher order terms in a and b are to be dropped, in accor- 
dance with the weak field approximation. In the above equations, k^ can 
be identified with the four-vector of the EM field before the arrival of the 
gravitational wave (that is, in flat space); its components are constant and 
satisfy the condition 

- 2k u k v + k 2 x + k 2 y = 0. (2.12) 

2.1 Short wave-length approximation 

In order to proceed further, set 



A 



a — CL a & , 



iS 



(2.13) 



where S is the eikonal function. Then Maxwell's equations F a J = imply 



+i 



+ a„ 



0. 



In the short wave-length approximation, the above equation is taken as an 
expansion in descending powers of S (see, e.g., Misner, Thorne and Wheeler 
[3J, Sect. 22.5). Defining the null four- vector 

K — R 

the first term in the expansion implies 

a a K a = (2.15) 
and the electromagnetic field turns out to be 

F aP = i{K a a p - K p a a )e lS . (2.16) 
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The second term in the expansion (I2.14p implies 

(a 2 K a ), a = 0, 
and the energy-momentum tensor, to first order, is 

T a /3 = a 2 K a Kp. 
In the present case, the eikonal function is 

S = + $, 

since 

crs,pS >v = o, 

again to first order in a and /3. Accordingly, 
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K + d u $ 

k v 

k x - k v (a"x + /3" y) - a'k x - p% 
\k y - k v {P"x - a"y) - p% + a'k y J 



(2.17) 

(2.18) 

(2.19) 
(2.20) 



(2.21) 



As for the four- vector a M , the components a x and a y must be independent 
of coordinates (though they may depend in general on fc M ), while a v = and 



ilt k v ((1 X K X + CLyKy 



(2.22) 



3 Stokes parameters 

It is convenient to define time-like and space-like unit four-vectors, t a and z a 
respectively, as 

t a = -L{-l + Itf, -1, 0,0} (3.23) 

Za =-L{-l-lH, 1, 0,0}, (3.24) 

such that t a t a = —1, z a z a = 1 and t a z a = 0. Then the frequency Q of the 
EM wave as measured by an observer at rest with unit velocity four-vector 
t a is 

n = -KX = -^[ R u + (1 + \h)K}. (3.25) 
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Clearly f2 depends on position and time. It is also convenient to define 

K\\ = K a z a , (3.26) 

and then ft 2 = K 2 + K 2 ± , where K ± = ^/K 2 + K 2 . 

Notice that in flat space, the frequency of the EM wave is just 

u = --j={k u + k v ) 

and the z component of the wave-vector is 

We can now define two four- vectors e« (i = 1,2) orthogonal to t a and 
K a : 

e« = KJHO, 0,-K y , K x } 

e (2) ° = -^z* + ^-{0,0, K t ,K y }. (3.27) 



The electric field in the frame at rest is defined quite generally as 

E a = F a %, (3.28) 

and accordingly the Stokes parameters can be constructed from the two scalar 
products E a 6a . Namely 

50 = \E a e^\ 2 + \E a e^\ 2 

5 1 = \E a e^\ 2 -\E a e^\ 2 

S 2 + iS 3 = 2(E a e^YE a e {2 \ (3.29) 

following the notation of Born and Wolf pi] . 

Now, 

E a = i[Vta a + (a^)K a }e iS (3.30) 

and since 

K a e {i)a = 0, 
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it follows that 

EJ® a = iVta a e® a e iS , (3.31) 

or explicitly 

E a e ma = i^{a y K x - a x K y )e lS , (3.32) 

E a e^ a = -i-^(a x K x + a y K y )e iS . (3.33) 

For algebraic calculations, it is convenient to express the Stokes parame- 
ters in matrix form as 

S + S, S 2 - iS 3 \ _ & kT f \a x \ 2 a x a* \ ^ 



S 2 + iS 3 Sq — Si J K 2 ' \ a* x a y \a y} 



where 



Ky K x 

-K x K y 



(3.35) 



The component K x and K y are given by (I2.21j) . 

The important point is that a x and (ly 3X6 constant, and therefore these 
quantities can be expressed in terms of their values before the arrival of the 
gravitational wave or, equivalently, in terms of the Stokes parameters s$ in 
flat space. The following relations are obtained with some straightforward 
algebra: 

q _ ^ 

tt 2 

S\ = — (si cos 25 — s 2 sin 25) 

n 2 

5 2 = —z(si sin 25 + s 2 cos 25) 

UJ 2 

53 = %s 3 , (3.36) 



CO 2 



where 



js 



-j—^(k x -ik y ){K x + iK y ). (3.37) 

The Stokes parameters can also be expressed in terms of two angles, % 
and ip, spanning the Poincare sphere (Born and Wolf [1]): 

51 = So cos 2\ cos 2tjj 

52 = So cos 2x sin 2ip 

5 3 = S sm2x- (3.38) 
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The angle ip specifies the orientation of the polarization ellipse and the angle 
X its ellipticity. 

The intensity of the EM wave before and during the passage of the grav- 
itational wave is given by the parameters sq and So respectively. From the 
definition of Q, it follows that 



n-u = — ^— {^[(fcX + M , 0^+(^/3 ,/ -A; y a , 0y] + a , (^-A;2) + 2^A; y ). 

(3.39) 

Similarly, it can be seen that the angle x remains constant, but the angle ip 
changes from ipQ to + ^ m the wake of the gravitational wave, where 



5 = {^[(^« ,/ -^/3'0x+(A; :E a /, +A;j / /3 /, )2/]-^(^-A:2)+2a / fc :c A; ?/ |. (3.40) 

It follows from the above formulas that the frequency of the EM wave 
is shifted from oj to Q. As for the polarization, a circularly polarized wave 
(si = S2 = 0) remains circularly polarized, but an elliptically or linearly 
polarized wave changes its direction of polarization. The above formulas 
show that the ellipticity of the EM wave is not altered, but the angle i/j, 
defining the axis of the ellipse, oscillates with the frequency of the functions 
a(u) and b{u). 



4 Summary of results 

We have obtained the relevant equations describing an EM wave in the pres- 
ence of a plane-fronted gravitational wave in terms of the eikonal function 

5 given by Eq. (I2.19P and two functions, a(u) and /3(u), characterizing the 
gravitational wave . 

The main result of the present article is that a gravitational wave inter- 
acting with a monochromatic EM wave rotates its axis of polarization by an 
angle given by Eq. (I3.40p . where the coefficients are the components of 
the four-vector defining the EM wave before the arrival of the gravitational 
wave; this wave four-vector changes to the form given by Eq. (I2.21[) in the 
presence of the latter. The rotation is oscillatory with the same frequency as 
that of the gravitational wave. 
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Appendix 

In this appendix, the EM field in the limit of flat space-time is worked out 
according to the formalism of the present paper. This permits to characterize 
the EM field before and during the passage of the gravitational wave and 
compare both cases. In flat space-time, the electric field is given by 

E = i(uti-a z k)e iS , (4.41) 

where k a = (u;,k) in Minkowski coordinates and 5" = k a x a . It is important 
to notice that the gauge used in this paper is such that a v = —a u = 0, and 
therefore a 1 = a z . Thus the condition k a a a = implies 

k x a x + k y a y = (lu - k z )a z , (4.42) 

and of course E • k = 0. 

The unit vectors 6a have purely space components: 

^ = kj_ ( k y ,k x ,0) 
e (2) = u-^^^ky-kj) (4.43) 

and thus 

. E = iukj 1 (k x a y — k y a x )e lS 
e< 2 ) -E = -iukj 1 (k x a x + k y a y )e iS . (AAA) 
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